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Abstract. We establish convergence in the diffusive limit from entropy 
weak solutions of the equations of compressible gas dynamics with fric- 
O | tion to the porous media equation away from vacuum. The result is 

D , based on a Lyapunov type of functional provided by a calculation of the 

00. relative entropy. The relative entropy method is also employed to es- 

tablish convergence from entropic weak solutions of viscoelasticity with 
memory to the system of viscoelasticity of the rate-type. 



1. Introduction 

The relative entropy method of Dafermos and DiPerna [3 El [10] provides 
an efficient mathematical tool for studying stability and limiting processes 
among thermomechanical theories. It is intimately connected to the second 
law of thermodynamics and has been tested in various situations involving 
stability and asymptotic behavior of shocks (e.g. [101 O H]), relaxation 
or kinetic limits in the hydrodynamic regime [271 Q] , stability and limiting 
processes among thermomechanical theories j5j EH H7J [5] . 

The method hinges on a direct calculation of the relative entropy between 
a dissipative solution and an entropy conservative (smooth) solution for the 
underlying thermomechanical process, which provides a remarkable stability 
formula [3 [6] . In more complicated situations involving the comparison of 
two solutions with shocks it is supplemented with additional information, 
e.g. [10ll3j[l9]. The objective of this article is to extend the relative entropy 
formula in situations where a dissipative solution of a thermomechanical 
system is directly compared to a dissipative solution of a limiting system. We 
use as test cases various paradigms of diffusive limits, the most significant 
perhaps being the validation of the limit from the Euler equations with 
friction to the porous media equation in the zero-relaxation limit. 

We consider the system of isentropic gas dynamics with friction 

1 

p t + - a\v x m = 

i <> ? i (1-1) 

m t + - div^ h -VxP(p) = — ,m 

e p e e A 

with the so called diffusive scaling, which captures the effective long-time 
response. In the limit e — > this system approaches the porous media 
equation 

Pt - A x p(p) = 0. (1.2) 



2 
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This problem has served as a paradigm for the theory of diffusive relaxation 
\24\ \18\ [12] and has been justified either by asymptotic in time analysis 
[IB"! [25l Ell El E], or via direct analysis of the relaxation limit, for weak 
solutions in [231 \22\ [24] or for smooth solutions near equilibrium in [H [20] . 

In this paper we compare directly a weak entropy solution of (II. lj) to a 
smooth solution of (|1.2j) using a relative entropy analysis (Proposition 12. ip . 
This, in turn, provides a convergence result to solutions of the porous media 
equation that stay away from vacuum (Theorems 12.71 and 12.8ft . The novelty 
of the present work is the simplicity of the proof following a Lyapunov 
type of analysis; in addition some new situations are analyzed (for instance, 
solutions approaching different end-states at ±oo), plus a rate of convergence 
is obtained. Finally, in the spirit of [21 [8], the relative entropy inequality is 
extended between entropy measure-valued solutions of the Euler equation 
and the porous media in Section [2. 4i 

We then test some other cases of diffusive relaxation using the relative 
entropy method. In Section [31 we consider the p-system with damping in 
Lagrangian coordinates and establish convergence to a parabolic equation, 
in the high- friction limit (Theorem 13. 3[) . In Section [H we consider the lim- 
iting process from viscoelasticity of the memory type (j4. 1[) to the system of 
viscoelasticity of the rate-type fj4.2|) in the diffusive regime. We provide a 
relative entropy estimation between the two theories and a convergence re- 
sult (see Proposition [4J] and Theorem l4.2p thereby extending for quasilinear 
systems previous convergence results in the semilinear case from [91 [11] . 

It is remarkable that in all those examples the dissipation of the approx- 
imating system can be split in two separate parts: the dissipation of the 
limit diffusion equation, and a second part that captures the dissipation of 
the approximating system relative to its diffusive-scale limit. 

2. ISENTROPIC GAS DYNAMICS IN EULERIAN COORDINATES WITH 

DAMPING 

We consider the system of isentropic gas dynamics in three space dimen- 
sions with a damping term: 

pt H — div x m = 

(2.1) 

1 m®m 1 . 1 

m t + - div x h -y x p{P) = — s-m, 

£ p £ £ z 

where t £ 1, i € I 3 , the density p > and the momentum flux m £ M 3 . 
The pressure p{p) satisfies p'(p) > which makes the system hyperbolic. 
An important particular case is that of the 7-law: p(p) = kp 1 with 7 > 1 
and k > 0. In (|2.ip . the variables (x,t) are already scaled in the so called 
diffusive scaling. In the diffusive relaxation limit e — > 0, solutions of (|2.ip 
formally converge to the porous media equation 

p t - A x p(p) = 0. (2.2) 
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The goal of this work is to study this limit via the relative entropy method. 

We recall that (rj,qi,q2,q3)(p,m) : M + x E 3 — > R x M 3 is an entropy- 
entropy flux pair for the hyperbolic system (12. ip if it satisfies the differential 
relations: 

dqj tdr\ 1 drj \ 1 dr\ 

= \7T + ~ m k-^ + - m j^\ — 

orrii \op p orrik / p orrii 



dqj 1 dr] ^7 



(2.3) 



T rrnmj- Vp'{p)- 



where i, j = 1, 2, 3, 5^ stands for the Kronecker symbol, and the summation 
convention is used. Moreover, the entropy r/(p, m) is dissipative (for the 
underlying relaxation process) if 

V( p , m )r/(p, m) • (0, -m) = - V m ry(p, m) • m < 0. 

An example of an entropy pair is provided by the mechanical energy 

1 1 777. 1 ^ 

r)(p,m) = -J—L + (2.4) 
and the associated flux of mechanical work 

1 |?7T/P / 

q(p,m) = -m — ^- + mh(p). (2.5) 

2 /T 

Here, /i(p) = pe(p), where e{p) is the internal energy of the gas connected 
to the pressure via e'(p) = ^*r. Accordingly, 

h »( p) = PM. p h '(p)=p(p) + h(p). (2.6) 

For the particular case of 7-law gases, h takes the form 

k 1 

-p 1 = rP(p) for 7 > 1; 



Hp) = { 7 - 1 7-1 

/cplog p for 7 = 1. 



Smooth solutions of (|2.ip satisfy the identity 

r/(p, m) t + - div x q(p, m) = — \v m r](p, m) ■ m = — < 0, (2.7) 

e e z e z p 

The mechanical energy rj(p, m) is dissipative for the relaxation process (|2.ip . 

2.1. Hilbert expansion. We start by reviewing the Hilbert expansion as- 
sociated to the relaxation process from (|2.ip to ()2.2j) . We introduce the 
asymptotic expansions 

p = p + ep\+ £ 2 P2 + • • • , 

m = mo + em-i + s 2 m2 + • • • , 

to the balance of mass and momentum equations in 112. lj) . and collect to- 
gether the terms of the same order, to obtain, respectively, 

0(e _1 ) div x m = 0, 
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O(l) d t po + div x m 1 = 0, 

O(e) d t pi + div x m 2 = 0, 



and 

0(e~ 2 ) m = 0, 

O(e^) -mi = V x p(p ), 

0(1) - m 2 = V x (p' (po)pi) , 

0(e) -m 3 = d t m x + div x ( mi (g3rril \ + \/ x (p'(p )p 2 + \v"(pa)p\ 

In particular, we recover the equilibrium relation mo = for the state vari- 
ables, Darcy's law ra\ = — V x p(po), and observe that po satisfies (|2.2p . 

Next, we focus on the asymptotic expansion of the entropy equation (|2.7p . 
and in particular on how the hyperbolic entropy (the mechanical energy) 
captures in the e — > limit the entropy structure of the porous media 
equation. Introducing the Hilbert expansion into (|2.7p and using ttiq = 0, 
we see that 

d t (h(p ) + eh'(po)pi) + - div x . (emih'(p Q ) + e 2 (m 2 h'(p ) + m 1 h"(p )pi)) 

\ m l\ 2 . f, , 2 Pl mi ' m A , rM 2\ 
= Ye mi ^ - 2 + 0(e ). 

Po V Po Po J 

Again, collecting together terms of the same order gives 

12 



0(1) h(p )t + div x (toi/i'(po)) = , 

Po 

0(e) d t (h'(p )p 1 ) + diyj; (m 2 h'(p ) + mih"(p )pi) = |mi| 2 ^ - 2^—^ 

Po Po 

Since toi = — V x p(po), the leading order term po in the diffusive limit sat- 
isfies the energy identity 

h(p) t - div x (h'(p)V xP (p)) = - |V ^ (P)|2 . (2.8) 



Equation (12. 8p captures the entropy dissipation of the porous medium 
equation (12. 2p and h(p) is indeed the entropy selected by Otto [26] in his 
gradient flow interpretation of (12. 2|) . 



2.2. Relative entropy identity. Let (p £ ,m £ ) be a weak solution of (12. ip 
that satisfies the weak form of the entropy inequality 

1 1 ItoP 
ij(p, m) t + - div x q(p, to) + < 0. (2.9) 

e e z p 

(We drop the e-dependence of (p e ,m £ ) except where emphasis makes it 
necessary.) Let p be a smooth solution of the porous media equation (|2.2j) : 
such a solution will also satisfy the entropy identity (|2.8p . Our goal is to 
devise an identity that monitors the distance between p e and p. 
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Such identities have been obtained via the relative entropy method in the 
context of problems of hyperbolic relaxation [17 1 127 } [l] . The relative entropy 
is defined as the quadratic part of the Taylor series expansion between two 
solutions (p,m) and (p, m); it takes the form 



rj(p, m\p,m) := r/(p, m) — ij(p, to) 
- V m r/(p, to) 



rj p (p,m)(p- p) 
rh) 



to 



1 



TP 



P 



m 
P 



+ h(p\p), 



(2.10) 



rrii 



while the corresponding relative entropy-flux reads 

qi(p,m\p,m) := qi(p,m) - qi(p,fh) - rj p (p,fh)(mi 

- V m ry(p, rh) • (/j(p, m) - fi(p, fn)) 

p p \ p p J p 

1, 2, 3, fi stands for the (vector) of the flux in (I2.ip 



1 

-TO; 

2 



where i 



fi(p,m) 



m \-p{p)h 

P 



(2.11) 



(2.12) 



and i, is the i-th column of the 3x3 identity matrix. 

Now the question arises about how to select to in (|2.10p . (|2.1ip . This 
relates to a significant difference among the hyperbolic relaxation and the 
diffusive relaxation frameworks: in the existing studies of hyperbolic relax- 
ation limits one compares an energy dissipative with an energy conservative 
solution. The fact that the limiting solution is energy conservative (and 
smooth) is an important restriction in the derivation of the relative entropy 
identities available in the hyperbolic relaxation framework (see [27} [I]). By 
contrast, by the nature of the diffusive relaxation framework, the solutions 
to be compared have both to be energy dissipative. To effect the comparison 
we select an e-dependent solution (p, to) that adapts itself in the relaxation 
process. 

A suitable selection of rh is proposed by rewriting (|2.2h in the form, 



- 1 i ■ - 

p t + - dw x m 

£ 







to 



-eV x p{p) 



(2.13) 



of the conservation of mass equation in (12.ip together with (a rescaled form 
of) Darcy's law. The energy identity (|2.8p may also be expressed in terms 
of (p, rh) as 



1 



h(p)t + - diva; (mh'(p)) 



1 \m\ 



e ' p 

In turn, (|2.13p is embedded into the system of Euler equations with relax- 
ation, plus additional terms purported to be higher-order errors. A simple 
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calculation shows that (p, m) satisfies 



Pt + ^d Xi mi = 
fh t + \d x Ji(p, fh) = --^m + e(p, fh) , 



(2.14) 



where (we use the convention of summation over repeated indices and) e is 



given by 

e := e(p, fh) 



1 I m®m 

- div x . = — 

e V P 



ed t V x p(p) 



^xP(p) ® V x p(p) 



EV x (p'(p)A xP (p)) 



(2.15) 



= e div x 

= 0(e) 
and is thus an error term. 

The main result of this section is: 

Proposition 2.1. Let (p,m) be a weak entropy solution of (|2.ip satisfying 
(|2.9p and let (p, fh) be a smooth solution of f|2. 13|) . Then, 



i](p,m \p,m)t + - div x q(p,m\p,m) < -^R(p,m\p,m) - Q - E, (2.16) 

e e z 

where 



R(p, m \p, m 
1 



m m 
P P 



Q 



-v 



n(p, m) 



Pxi 



E = e(p, m 
and e(p, fh) is defined in (|2.15p 



p I m m 
P\P P 





fi(p,m\p, fh] 



(2.17) 



Remark 2.2. The following remarks are in order, concerning the terms ap- 
pearing on the right of ()2.16|) : 

(a) The coefficient of the quadratic term Q depends only on (p,fh); it is 
explicitly given by 



1 



(p, m)p Xi + rj mkrrij (p, m)d Xi m k 



1 



rn , 



-d XiXj h'{p) , 



and is thus of O(l) in e. 

(b) Since e(p,fh) = 0(e), the coefficient of the error term E is of 0(e). 

(c) The term R(p, m\p,fh) captures the dissipation of the relaxation system 
(12. ip relative to its diffusive scale limit (12. 2p . It turns out to be the quadratic 



RELATIVE ENTROPY IN DIFFUSIVE RELAXATION 



part of the dissipative relaxation term with respect to (p,m). Indeed, for 
R(p, m) = V m r] ■ m = R p = rj pmj mj, R mi = r) mi + T} mim .fhj. Using 

1 



m 



rrin 



V m rj(p,m) = — , n pmj {p,m) = f, r) mkm Ap,m) = -5 kj , 




P 

we compute the Hessian of R, 
and see that it has eigenvalues 

2 2 2 

Ai = 0, A 2 ,3 = - > 0, A 4 = - + ^tH 2 > > 
P P P 6 

and is thus positive semidefinite for p, p > 0. A direct computation also 

2 

and justifies the notation used in 



p 4 

mi _ . 

-2-^- -o tJ 

P 2 P 



*J=1,2,3 



m m 

P P 



shows that R(p,m\p,m) = 

(|2,17P i . The property R(p,m \p, fh) > is instrumental in the forthcoming 
stability analysis; a similar property holds for all examples treated in this 
article. 

Proof of Proposition \2.1\ By hypothesis (p, m) satisfies the weak from of the 
entropy inequality 



11 1 

n(p, m) t + - div x q(p, m) < ^V m 7?(p, m) ■ m = — 5- 

e e z e z p 

Also, (p, fh) satisfies the energy identity 



(2.18) 



r)(p,m) t + - div x q(p,m) 



1 



VrnViP, m) -m + V m 7](p, m) ■ e(p, m) 



1 Iml 



e 2 p 



+ V m rj(p,m) ■ e. 



From ([2~T|) and (|2~Tlj) we obtain 



(p - p) t + ^d Xi (rrii - mi) = 

(m - m) t + ^d Xi (fi(p,m) - fi(p,fh)) = -^(m - fh) - e 
and use the smoothness of (p, fh) and (|2,14p to compute 



(2.19) 



(2.20) 



dt[r] p (p,m)(p - p) + V m n(p,m) ■ (m - m) 



+ -d Xl [rip(p,m)(mi - rrn) + V m rj(p,m) ■ (fi(p,m) - fi(p,m] 



-—~{m — m) — e 



V m r)(p,m) 



+ 9 t (v P {p, m)) (p- p) + d t (V m r/(p, m)) • (m - m) 
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+ -d Xi (rj p (p, m)) (mj - urn) + -d Xi (V m r](p, m)) ■ (fi(p, m) - fi(p, m) 
-~2 V m ry(p, fa) ■ (m-fa) - V m r)(p, fa) ■ e 



m 



Pt 



p-p 
m — fa 



P_Xi 



rrii — rrii 
fi(p,m) - fi{p,fh) 



--2^mV(P, m) ■ (m-m) - V m r](p, ra) ■ e(p, ra) 



+ V(p,m)V(P,rn) (_l 
1 







+ -V 



(p,m 



)V(P, 



m 



■^ra + e(p,m) 

Pxi 



p-p 
m — fa 







fi(p,ra\p,m) 



(2.21) 



=: J. 



To obtain (j2.21[) we used the identities 



^t P ,m)V(P^) 



P-P 
m — fa 



V 



(p,m) 



rj(p, ra 



V (p, m )V(p,m) 
V 2 M r]{p,rh) 



d Xi rrii 
d x Ji(p,m 

5n 5 i2 S i3 

V(p,m)fi(p,rn) 



Pxt 



m 



Xi 



in 



Pxi 



V( p , m )/i(p,m) 

rai — rhi 
V( p , m )/i(p,m) • (m -fh) 



P-P 
m — fa 

P-P 
ra — fa 



(resulting from the entropy consistency relations) and the notation 
fi{p, m\p, fa) = fi(p, ra) - f { (p, fa) - V( p>m )fi(p, fa) ■ {p - p,m - fa) . 
Finally, combining ([2TTU]) . ([2TTD . ([2TT8]) . ([2TT9]) and ff272T^ . we conclude 

ri(p,m\p,m)t + - div x q(p,m\p, fa) 



< 



V m r?(p, ra) - ra - VmViP, m) ■ m - \7 m 7](p, m) ■ (m - m) 
- r] pmj (p,fn)fhj(p - p) - ij mkmi (p,fh)fhi(rn k - fa k ) 



-V 



(p,m 



)V(p,m) 



Pxi 

m Xi 





fi(p,ra\p,m) 



ei(p,ra)[rj pmi (p,ra)(p - p) + v. mkmi {p,ra)(ra k - m k ) 



--: -^D-Q-E, 
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where 



/ mi , , 1. . \ , _ . p I m m 

E = ei[p,m) — -y(p - p) + -(mi - mi) = e(p,m) ■ - 

P z P ) P\PP 

-D = - ■ (m-m) + -zs-(p - p) —m e ■ (m k - m k ) 

P P P 



12 

m\ 


1 - 12 

\m\ 


P 




P 


m 


m 


2 


P 


P 





= p 

and Q is as in ([2717]) . □ 

We conclude this section with two lemmas. The first establishes (under 
additional hypotheses on p) a bound of the quadratic term Q in terms of 
the relative entropy (|2,10p . 

Lemma 2.3. Assume that p(p) satisfies for some A > 

p"(p)<A^- Vp>0. (A) 

Then h(p) in (|2.6p verifies, for some c > 0, 

p(p\p)<ch(p\p) Vp,p>0. (2.22) 
f/iere exists a C > suc/z f/mf /or any i = 1,2, 3, 

\fi(p,m\p,m)\ < Cr,(p,m\p,m). (2.23) 

Proof. Recall that h" = Hypothesis ((23) implies p" < Ah" . One easily 
checks the identity 

P(p\p) = P(p) ~ P(p) ~ p'(p)(p ~ P) 

= {P-P? I I P"(sp + (1 - s)p)dsdr , 



JO JO 

and a similar identity holds for h(p\p). Then (|2.22p follows from p" < Ah" . 
A direct computation using (|2.12p shows that, for i, j = 1,2,3, 

rmmj fhifhj fhifhj , 

fij(p,m\p,m) = ^ + —^(p- p) 

P P P 2 

- -zi^uffij +fhiSje)(me - m e )) +p{p\p)5ij 



/ in, n H > ; in, ill-, > . . _. 

= P — - -rr )+P(P \p)Sij , 

VP P J V P P / 

and ([2722]) gives ([2723]) . □ 

The second lemma indicates a relation between the "metric" induced by 
the relative entropy (I2.10P and more traditional norms. 
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Lemma 2.4. Let h G C°[0,+oo) n C 2 (0,+cx)) satisfies h"(p) > and 

h(p) = — k —p< + o{p>) , asp^+oo (2.24) 
7-1 

for some constant k > and for 7 > 1 . If p G K = [5, M] with 5 > and 
M < +00, f/ien there exist positive constants Ro (depending on K) and C\, 
C2 ( depending on K and Rq ) such that 



Kp\p) > 




p\ 2 , /or < p<R , p G K, 
pp, for p> Rq, p G K. 



Proof. Since p G K, there exist positive constants A and B such that for 
any < p < +00 

h(p \p) = h(p) - h(p) - h'(p)(p -p)> h(p) -A- Bp, 

Hp\p) > Hp) -a- Bp 

\p-p[r- P~< 
In view of (|2.24p . there exists i?o depending on K such that 

1 k 

hip Ip ) > -\p — pP , for p > Ro and p G K. 

2 7 — 1 

Consider next the positive function 

B(p,p):=P^L, P e[0,Ro],peK,p^p 

\p-pr 

and note that lim p ^p B(p, p) = > and thus B(p, p) G C°([0, i? ] x A"). 
Hence, we have 

h(p \p) > m|p — /?| 2 for /) G [0, i2o] and p £ K 

where m := min pe[0;il . ] ; p eA ' B(p, p) > 0. □ 

Remark 2.5. Hypothesis (jXJ) holds for large classes of pressure laws. Indeed, 
is trivially satisfied for concave increasing pressures. For convex pres- 
sures, one checks that it holds for p(p) = kp 1 with 7 > 1. By contrast, (JSJ) 
is violated for p(p) = e p . For a 7-law gas with 7 > 1, one computes that 
h(p) = ^tP(p) an d thus (I2.22P holds as an equality. 

Regarding the growth condition (I2.24L if p G C°[0,+oo) n C 2 (0,+oo) 
satisfies p'(p) > and 

p'(p) = k 1P ^- x + oip 1 - 1 ) , as p ->• +00 , (B) 

with fc > 0, 7 > 1, then h{p) defined through h"{p) = verifies hypothesis 
(|2.24p and the results of Lemma 12.41 apply in that case. 
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2.3. Convergence in the diffusive relaxation limit. Proposition 12. II is 
used in order to prove convergence from the Euler equations with friction 
in the diffusive limit towards the porous media equation. We carry out the 
analysis in two frameworks: 

• multi-d periodic solutions; 

• 1-d solutions in the real line with (possibly distinct) constant states 
p± at ±oo. 

In both cases, the main hypothesis is that p is a smooth solution of (|2.2p 
that sits away from vacuum. 

Remark 2.6. It is worth to observe that other possible frameworks can be an- 
alyzed with these techniques, and we restrict our ourselves to the aforemen- 
tioned cases to avoid further technicalities. For instance, with small mod- 
ification in the arguments below, we can consider multi-d solutions (p, fh) 
such that p — > > as \x\ — > +oo; m = —eVp{p) and such that p > 0, 
p-p*, m e L\[0,T] x R 3 ). 

2.3.1. Multidimensional periodic solutions. In the periodic case, we 
work within the following framework, collectively referred to as (Hi): 

(i) (p , m) : (0, T) x T 3 — > M 4 is a (periodic) dissipative weak solution of 
(|2.ip with p > 0, satisfying the weak form of (12. 1|) and the integrated 
form of the entropy inequality (j2.9|) : 

|2 \ i l™|2 



[0,+oo)xT 3 



(1 \tt) I \- 1 Im I 

-— + h(p) 6{t) - -z^0{t) dxdt 
2 p ) e 2 p 

+ / (I^T + Hp)) t A0)dx>0, (2.25) 
y T 3 \2 p j t=o 

where 9{t) is a nonnegative Lipschitz test function compactly sup- 
ported in [0, T). The family (p £ ,m £ ) is assumed to satisfy the uni- 
form bounds 



sup / p e dx < Ki < oo , (2.26) 

£€(0,T) JT 3 

sup / h h(p E ) dx < K2 < 00 , 

te(o,T) Jt 3 2 P £ 

which are natural within the given framework, and follow from cor- 
responding uniform bounds on the initial data, 
(ii) p is a smooth (C 3 ) periodic solution of the multidimensional porous 
media equation (12.2j) that avoids vacuum, p > c > 0; in is defined 
via fh = —eVp(p), 
The function 

fit) = / r]{p } m\p,fh)dx (2.27) 



JT 3 

will be used as a measure to control the distance between two solutions. We 
prove: 
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Theorem 2.7. Let T > be fixed and assume p(p) satisfies (|A|) and (iBj) , 

Under hypothesis (Hi), the stability estimate 

<p(t) < C(v?(0) + £ 4 ), te[0,T], (2.28) 

holds, where C is a positive constant depending only on T, K\, p and its 
derivatives. Moreover, if ip(0) — > as e 10, then 



sup ip{t) —7-0, as e 10. 
te[o,T\ 



(2.29) 



Proof. We proceed to establish the integrated version of (|2.16p under the 
regularity framework (Hi). To this end, we introduce in f)2.25j) the choice 
of test function 



6{r) := { 



Taking the limit k \, gives 



1, for < r < t, 

*=T + 1, fort<T<t + K, 
0, for r > t + K. 



(2.30) 



1 m 



2 p 



+ h(p) I dx 



< 



1 



r=0 



m 



-dxdr . 



£ Jo Jts P 



Next, integrating (|2.19|) over (0, t) x T 3 , gives 



1 \m\ 



3 \2 p 



+ h(p) I dx 



f f f 1 Imp m , , , 
< / / — ^ — = — h — ■ e ) axar . 



Finally, to justify the calculations leading to (|2.2ip . we start from the 
weak form of (I2T20]) : 



0t(p - P) + -4>xi (m - mi) dxdt 



(x,0)(p-p) 



t=o 



dx = , 



(2.31) 



if) t ■ (m - m) + -ip Xi ■ (fi(p, m) - fi(p, m)) dxdt 



ip(x, 0) • (m — m) 



dx 



t=o 



ra — m) — e I dxdt , (2.32) 



where (j), ip are Lipschitz test functions compactly supported in [0, T) x T 3 
and ip is vector valued. Using the test functions 

4> = 9(r)uj(x)r}p(p, fh) , V = 9{t)u{x) V m n(p, m) 

with 9{t) as in (I2.30P and w(x) = 1, and upon taking k J, 0, this gives 



T 3 



n p (p,m)(p - p) +V m r](p,m) ■ (m - m))dx 



r=0 J JT3 



J dxdt , 



(2.33) 
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where J is as in (I2.21|) . Combining the above inequalities leads to 



e Jo JY-s 



R(p, m \p,m)dxdT < (p(Q) 



o Jt 3 



(\Q\ + \E\)dxdT (2.34) 



where Q, E and R are given in (|2.17|) . 

Using ([237|> . Remark O (a), Lemma and (l2T2Tjl . we deduce 



o JT 3 



where Cj depends on llc^x^j ^'(p) \\l°° • The error term E in (12. 17j) is esti- 
mated by 



\E\dxdr < — 
^T 3 2 



o Jt 3 



and, by ([21511 and (12251) . 



o .7t 3 



pdxdr + 



1 









m 


m 


/'/" 






/O JT 3 


P 


P 



pdxdr <C\£ t . 



dxdr 



(2.35) 



where C2 depends on K\ , T and p through the following norms of derivatives 
up to third order: 



— div x 
P 



V x p{p) g VsP(p) \ 
P / 



+ 



-d t V x p(p) 

p 



Introducing the above estimates into (|2.34p . we obtain 

<P(*) + I I R(p,m\p,rh)dxdT 
2e 2 Jo Jt 3 

< ip(0) +Ci f ^(r)dr + C 2 e 4 t • 
Jo 

Gronwall's inequality then implies 

<p(t) < C(<p(0) + e 4 ), t£(0,T] 

and (I2T29]) follows. 



□ 



2.3.2. The Cauchy problem on the real line. Next, we consider the 
Cauchy problem in one-space dimension for 



pt + ~m x = 

e 



m t + 



1 / 771 



1 



(2.36) 



-m . 



+ P(p), 
£ \ p J x E< 

To avoid unnecessary technicalities with the behavior as |x| —> 00, we assume 
the initial data (po,mo) take constant values outside a compact set [— Rq.Rq], 

(p (x),m (x)) = (p-,0) forx<-i? , 

(p (x),m (x)) = (p + , 0) for x > R 
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for some p± > 0. By the finite speed of propagation property, any solution 
(p, m) will assume the same values outside the cones x < —Rq — kt and for 
x > Rq + kt, respectively, with k calculated in terms of the maximum wave 
speed on the range of the data. 
Let p > be a smooth solution of 

Pt-P(f>)xx = (2.37) 

with initial data po taking constant values 

Po(x) = P- for x < —Rq , po{x) = p + for x > Rq , 

outside some compact set \—Rq,Rq\, with p± > as above. By standard 
theory for the porous media equation (see [28]), the solution of (|2.37|) satis- 
fies p(x,t) > c > 0, and satisfies p(x,t) — > p± as x — > ±oo with sufficiently 
fast decay (in fact exponential). Defining to = —ep(p) x , we obtain m — > 
as x — > ±oo. 

By modifying the entropy pair (|2,4p - (|2.5p (using a trivial linear pair), we 
define 

fj(p,m) = V (p,m) - h ^- h ^ (p-h p+ + p _)) - hh(p+) + h(p-)), 

p + — /0_ V I /I 

~t \ i \ h (p+) - Hp ) 

q{p,m) = q{p,m) m 

P+ ~ P- 

so that fj(p±,0) = 0. The resulting (fj — q) is an entropy pair, and vanishes 
at the end states (p±,0). 

We next summarize the framework (H2) for the relaxation limit: 

(i) (p,m) : (0, T) x R -> R 2 with p > is a dissipative weak solution of 
(|2.36p . that is, it satisfies the weak form of (I2.36P and the integrated 
form of the entropy inequality 

fj(p,m)6(t)dxdt + / 7/(/j, to) 6(0)dx 

[0,+oo)xR Jr t=0 

1 ff ItoI 2 
> — / / i—i-e(t)dxdt, 

£ ii[0,+oo)xK P 

with 9{t) a non negative Lipschitz test function compactly supported 
in [0, +00). The family {p £ ,m £ ) is assumed to satisfy the uniform 
bounds 

sup (f \p £ -p-\dx+[ \p £ - p + \dx] < Ki < 00 , (2.38) 
te(o,T) \J-00 Jo J 

sup / f](p £ , m £ ) dx < K2 < 00 , 
te(o,T) 

with K\, if 2 independent of e. Of course, this dictates analogous 
uniform bounds on the energy norm of the initial data (/9q,TOq). 

(ii) p is a smooth (C 3 ) solution of (|2.37[) that satisfies p > c > 0; to is 
defined via to = —eVp{p). 
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We now denote by 



fj(p, m \p,m)dx. 



(2.39) 



This will replace (|2.27[) as a yardstick for measuring distance between solu- 
tions in the one-dimensional Cauchy problem. Then we have: 



Theorem 2.8. Let T > be fixed and assume ([A]) and ((B|) hold. If (p,m) 
and (p,fh) are under the framework (H2), the following stability estimate 
holds: 

m<c(m+e A ), te(o,T\, 

where C is a constant depending only on T, p±, p and its derivatives up to 
third order. Moreover, if (j){0) — » as e -10, then 

sup 4>{t) — > 0, as e \,0. 
te[o,T] 

Proof. Proceeding along the lines of the proof of Theorem 12.71 one derives 
the analog of (|2.34p for <j> in (|2.39j) . There is however a difference in the 
derivation as applies to the Cauchy problem: the equations (|2.3ip and (|2.32p 
hold for test functions compactly supported in [0, T) x M. Thus we introduce 
the test functions 



tp = 9(T)uj(x)rip{p,m) . 
where 9{t) defined in (|2.30[) and 

' 1 for 



ip = 9{T)uj{x)r] m (p,m) 



Uj(x) 



1 + 
1 + 





R-x 
S 

R+x 



R < x < R, 
for R < x < R + 5 
for - R - 5 < x < -R 
for x > R + 5 orx < —R 



5, 



into (|2.3ip and (|2.32p . Sending R — > 00, using the asymptotic properties in 
x of (p, in) and (p, fh), and subsequently sending k J. 0, we obtain the analog 
of (12331) and through that the analog of fl!TMj) . 

A second difference lies in replacing (|2.35j) by the estimation 



pdxdr < 



(p-l x<0 + P+h<o) I dxdr 



+ max{p_,p + } 
< CeH , 



dxdr 



where we used (I2.15P and the constant C depends on T, K\ in (12.381) . and 
also on p through the L°° norms of space-time derivatives up to third order 
and the norms 



div-r 



p 



Lf (L2(R)) 



+ 



dt^xP(p) 



P 
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Again using Gronwall, we deduce 

Ht) <C(<X0) + e 4 ), te(p,T\ 

which completes the proof. □ 

2.4. Relative entropy for entropic measure- valued solutions. A vari- 
ant of the relative entropy identity can be derived for comparing entropic 
measure- valued solutions of (|2.1j) with smooth solutions of (|2.2p . Such cal- 
culations are in the spirit of the recent works [8] , the difference here being 
that two dissipative systems are compared. 

Let v = {v x ,t] r x t)^Q T ^ e a parametrized family of probability measures 
(Young measures) that acts on continuous functions f{X p ,X m ), (Ap,A m ) G 
R + x R 3 , via 

{vx,t,f) = J /(Ap, X m )dv(X) 

and such that the integral (when defined) is measurable in (x,t) € Qt- A 
measure- valued solution of (12. ip consists of a Young measure {y x j}, ^ g g T 
with averages 

Wx,t, Ap) = p , (v X)t , A m ) = m , (2.40) 
that satisfies in the sense of distributions the measure- valued version of (|2.ip 

1 

Pt + ~ diva; m = 
m t + - div x (u xt , ^ m f Am + p(Xp) I) = — \m . 

The Young-measure v = {yx,t\ r x t)eQ T * s canec ^ an en t ro P v measure valued 
solution if it also satisfies in the sense of distributions the averaged version 
of the entropy inequality 

1 1/1 Am 1 2 \ 

dt{v x ,t , f?(Ap, Am)) + - div^i/^ , q(X p , X m )) < —^(v X)t , — ^— J , (2.41) 
for i] - q as in ([2T4>([X5]) . 

Proposition 2.9. Let v = {v x ,t}t x t)eQ T satisfying (|2.40p be an entropy 
measure-valued solution of (|2.1|) . and let (p,fh) be a smooth solution of 
(|2.13p . Then, we have the following averaged relative entropy inequality 

dt{v x ,t , r)(X p , X m \p,m)) + - div x (v Xjt , q(X p , X m \p,m)) 

< £ 1 / . A TO fh 2 ^ ( 2 - 42 ) 



A P p 



where 



Q - -V 2 {p ,m)V(P,rn) • ((^,/ j( Ap°A m |p,m)>) ' 
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_ e(p,fh) I (X m fh 

£- — = — ' \ v x,ti A p — — 

P \ V x p P 

and e(p,fh) is defined in (|2.15l) . 

Proof. We use (|2.10p to define the averaged relative entropy 

(v x>t , r](X p , X m \p,fh)) = (u X)t , r)(X p , X m )) - rj(p,fh) 

- 7] p (p, fh){u x>t , X p - p) - V m 7](p, fh) • (u Xit , X m - fh) . 

The inequality (|2.42p is built by using (|2.4ip , (|2.19p and the averaged version 
of (|2.20p and following verbatim the steps and calculations in the proof of 
Proposition 12.11 □ 

3. The ^-system with damping 

The p-system with damping in one space dimension is the system of con- 
servation laws 

1 

u t v x = 

1 £ 1 ^ 

V t t(u) x = —rV , 

e e z 

where r satisfies t'[u) > to guarantee strict hyperbolicity. The system 
(|3.ip is a model either for elasticity with friction or for isentropic gas dy- 
namics in Lagrangian coordinates (denoted by (x,t)). Then u stands for 
the strain (or the specific volume for gases), v for the velocity and r for the 
stress. 

In the high friction limit e — > 0, solutions of (|3.ip converge towards a 
solution of the parabolic equation (see [23] ) 

u t -T(u) xx = 0. (3.2) 

We will indicate in this section a simple proof of that convergence using the 
relative entropy identity. 

For concreteness, we interpret (|3.ip as a model for shear motions, n, v 
take values in R. We place the hypothesis that r : R — > M satisfies t'(u) > 
and the growth assumptions 

t(u) = ±\u\ p + o(\u\ p ) , as n — > ±oo , (H) 

for some p > 1. 

3.1. Preliminaries. The approach uses the mechanical energy 

1 f u 

£(u, v) = -v 2 + W(u) , where W(u) = / T(s)ds 

2 Jo 

is the stored energy. The associated flux is 

T{u, v) = —vt(u) 
and they satisfy the entropy inequality 

£(u, v) t + -T(u, v) x <-\v 2 <0 (3.3) 
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indicating the dissipation of the mechanical energy. 

The minimum of the mechanical energy £ (u, v) on the "equilibrium man- 
ifold" of the relaxation process M. = {(u, v) : v = 0} is achieved and is given 
by 

W(u) = £(u,0) = mm£(u,v) . 

M 

Moreover, solutions of (j3.2[) satisfy the following energy estimate: 
or equivalently 



£(u,0) t = {t(u)t(u) x ) x - (t(u) x ) 2 , 



£{u,0) t + F(u,T{u) x ) x = -{t{u) x ) 2 . (3.4) 

Relation (13. 4h captures the equilibrium version of (I3.3H . as can be seen by 
applying the Hilbert expansion to the relaxation system (]3.ip . 
Indeed, introducing the Hilbert expansion 

u = uq + eu\ + e 2 U2 + ... 

^ = ^0 + 6^1+ e 2 v 2 + ■ ■ ■ 
to (13. ID . we see after collecting the terms of similar orders that 
O^ 1 ) d x v = 0, 

0(e°) d t u + d x vi = 0, 

and 

0(e- 2 ) v Q = 0, 

0(e _1 ) vi = p(u ) x , 

O(e ) v 2 = (p'(u ) Ul ) x . 

In particular, we recover the equilibrium relation vq = 0, the Darcy's law 
v% = t(uo) x , and the diffusion equation (|3.2p satisfied by uq at equilibrium. 
If the same expansion is introduced in (|3.3p . we obtain 

£{u ,0)t - (t(uq)vi) x = -v\ 

which yields (I3.4p upon using Darcy's law v\ = t(uq) x . 

3.2. Relative entropy estimate and study of the relaxation limit. 

To analyze the relaxation process, we consider the quadratic part of £{u,v) 
with respect to the "algebraic-differential equilibrium" (u, v), where u = uq 
and v = ev\ = et{u) x . Namely, 

£(u, v\u,v) = £(u, v) — £(u, v) — £ u (u, v)(u — u) — £ v (u, v)(v — v) 

= 7^( v ~ v) 2 + W(u \u) . 
As corresponding flux we shall consider 

T(u, v\u,v) = T(u, v) — T(u, v) + £ u (u, v)(v — v) + £ v (u, v)(t(u) — t(u)) 
= —{v — v)(t(u) — t(u)) . 
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As in the previous section, to simplify the calculations, we rewrite the equi- 
librium equation (13. 2\ as follows: 

( _ 1_ 

u t v x = 

(3.5) 

- 1 / X 1 
V t T(U) X = ttV + V t . 

£ £ A 

In this way, we are able to treat the term vt = £r(u) x t as an error of order 
0(e). A direct computation, along the lines of Proposition [2~T1 gives: 

Proposition 3.1. For any weak, entropy solution (u,v) of (|3.ip and any 

smooth solution (u, v) of (]3.5p it holds: 

£(u,v\u,v) t + -J r (u,v \u,v) x < — \{v-v) 2 +t(u) xx t(u \u)-£T(u) xt (v-v) . 

£ £ A 

(3.6) 



The terms in the right hand side of (13.6P are analogous to the terms in 
(|2,16j) of Proposition 12.11 for the Eulerian case, namely, the first term is 
dissipative and is due to the damping of the relaxation system relative to 
its diffusion limit, the second is quadratic in the flux, and the last term is 
a linear error term. The quadratic term is estimated with the help of the 
following lemma from [8]. 

Lemma 3.2. Let r G C 2 (M) satisfy t'(u) > and (|H|) . If u takes values in 
a compact set K, there exists constant C such that 

t(u\u) < CW(u\u) , VueR, ueK. 

Proof. Since u G K, K a compact set, using (jH]) there exists a constant C 
such that 

\r(u\u)\ < C(\u\ p + 1) , u E K . 



Moreover, 



WM = Um r(«) _ i 



|u|-H-oo \u\^+oo (p + l)|n| p sgn(n) p+1 

and therefore, for some c and A, we obtain 

W{u \u) = W(u) - W(u) - W\u)(u -u)> c\u\ p+l - A , 

for u G R, u € K. We select U$ so that for \u\ > Uq we get c|u| p+1 — A > 
C(\u\ p + 1). Then 

\t(u\u)\<CW{u\u) V\u\>U ,ueK. 

On the complementary interval u € [— t/o, J7o]> we have 

t'{u) > c> 

and thus 

t(u\u) <Ci(u- u) 2 < C 2 W(u\u) V |u| < U , u£K. 

□ 
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Using Proposition l3.14 we obtain the main stability and convergence result 
in terms of the quantity 

$(t) = / £(u,v\u,v)dx. 
Jr 

Theorem 3.3. Assume r satisfies t' > and dHJ). Let u be a smooth 
solution of (|3,2p . defined on Qt = Mx [0, T) with T > fixed, v = er{u) x , 
be such that u is bounded and \\t(u) x i)\ L 2iq t \ < K < oo. Let (u,v) be a 
weak, entropy solution of (|3,ip such that 3>(0) < +oo and 

T(u, v \u, v) — > 0, as x — > ±oo. 

T/ien i/ie following stability estimate holds: 

$(t) < C($(0) +e 4 ), tG[0,T), (3.7) 

where C is a constant depending on T, the properties of t(u), and the func- 
tion u and its derivatives. Moreover, z/<3?(0) — > as e J, 0, then 

sup ->■ 0, as e 10. 

te[o,T] 

Proof. The proof proceeds along the lines of Theorem 12.71 and Theorem 12.81 
here we shall just sketch it. 

We integrate ([M]) over R X [0,t], t < T. The right hand side of is 
estimated using Lemma 13.21 

\t(u) xx t(u \u)\ < C\W{u\u) < C£(u,v \u,v) , 

and Young's inequality 

1 £^ 

\eT(u) xt (v - v)\ < -^{v - vf + —\r{u) xt \ 2 ■ 
Then, (|3,7p is a direct consequence of the Gronwall Lemma. □ 

4. VlSCOELASTICITY WITH MEMORY 

We conclude with an example where the diffusive scaling limit is a hy- 
perbolic - parabolic system. Consider the following 3x3, one dimensional, 
quasilinear system of viscoelasticity with memory effects: 

u t - v x = 

/ \ 1 

v t - a(u) x - -z x = ( 4 _-q 
fi 1 

Z t V X = yZ, 

where /i > and the elastic stress function a satisfies the usual condition 
o~'{u) > which guarantees hyperbolicity. The above system describes a 
one dimensional viscoelastic material for which the stress S = o~(u) + is 
the sum of an elastic part and a viscoelastic part of the memory type (see 
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(|4.ip <t). The system is scaled appropriately so that it relaxes as e — > to 
the equations of viscoelasticity of the rate type, 

ut - v x = 
v t - a(u) x = fiv xx . 

In the latter system, the total stress T = a(u) + fiv x consists of an elastic 
part and a Newtonian viscous stress. We refer to [9J |TT] for studies of a 
corresponding semilinear relaxation framework, using energy bounds. Here, 
we focus at the quasilinear level, and pursue a relative entropy analysis to 
explore the relation between the two systems. 
The mechanical energy for (|4.ip is 

f u 11 11 

E(u,v,z) = J a{s)ds + -v 2 + —z 2 = ^{u) + -v 2 + —z 2 , 

with energy flux 

¥ £ (u, v, z) = —(ea(u)v + vz) . 
Weak solutions of (14. lh are required to satisfy the entropy inequality 

E(u, v, z) t + -F e (it, v, z) x < — \z 2 , (4.3) 

manifesting the dissipation of the mechanical energy. (Smooth solutions of 
(|4.ip satisfy (|4.3p as an equality.) 

The equation capturing the dissipation of mechanical energy for smooth 
processes of viscoelasticity of the rate type (|4.2p reads 



E(«, v, 0) t + Fi(«, v, <t{u) x ) x = -fi{v x ) 2 , (4.4) 

where 

E(u,v,0) = E(w) + ^v 2 

is the equilibrium energy for E(u, u, ^) and 

Fi(«,u,<T(it)a.) = (cr(u)f + fj,vv x ^). 

Note that (j4.4j) is the leading order (with respect to the relaxation parame- 
ter) asymptotic development of the energy dissipation inequality (|4.3[) . This 
may be seen, as in the previous sections, by expanding (|4.3p in terms of the 
Hilbert expansion; we omit the details here. 



4.1. Relative entropy estimate and study of the relaxation limit. 

Following the general procedure, outlined in Section \2.2\ we recast the equi- 
librium system (|4.2p and the corresponding stress-strain response in the 
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variables (u, v, z) with z = eliv x as follows: 

u t - v x = 

1_ 

v t - cf(u) x z x = 

e 



(4.5) 



A* _ 
^ 

e 



1 _ 



where we shall treat the term zt as an O(e) error: 

= e/iw xi = epb{a(u) x + AwW)^ . 
We define the relative entropy and relative entropy flux, respectively, 

E(ii, v, z\u,v,z) = E(u, u, z) — E(u, t>, z) 

— E u (u, u, z)(u — u) — ¥. v (u, v, z)(v — v) — ¥ z (u, v, z)(z — z) , 
¥ e (u, v, z\u,v,z) = ¥ £ (u, v, z) - ¥ e (u, v, z) - ¥ u (u, v,z)(- e(v - v)) 

— E v (u, v,z)( — e(a(u) — a(u)) — (z — z)) — E z (u, v, z) (y — v) , 
and derive the relative entropy identity: 

Proposition 4.1. Let (u,v,z) be a weak entropy solution of (14. f p and let 

(u,v,z) be a smooth solution of (14.50 . Then 



dt¥j(u, v, z \u, v, z) H — d x ¥ £ (u, v, z \u, v, z) 

< n (z - z) 2 + v x a(u \u) - ev xt (z - z) . 

Proof. We use (|4.3p and rewrite (|4.4p as follows: 

1 /- -X 1 

-¥ e (u,v,z) x = 

e Lie 

Moreover, a direct computation shows 



(4.6) 



E(u, v, z) t + ^F e (u, v, z) x = -^z 2 + ezv xt . 



dt (cr(u)(u — u) + v(v — v) H — z(z — z)^j 

H — d x ^ — ea(u)(v — v) — ev(a(u) — cr(u)) — v(z — z) — z(v — v)) 

1 __ __ _1_ 

= ~z(z - z) - ezv xt + a {u)ut(u - u) + v t (v - v) + —z t (z - z) 

fie^ ll 

1 _ _ 1 

- a(u) x {v - v) - v x (a(u) - a(u)) - -v x (z - z) - -z x (v - v) 
2 

= t;z{z — z) + ev x t(z — 2z) — v x (a(u) — a(u) — a'(u)(u — u)) . 

Finally, putting all relations together we obtain (14. 6p . 



□ 
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Proposition 14.11 suggests to measure the distance between systems (|4.1|) 
and (|4.2|) via the quantity 

Vl/(i) = / E(it, v, z \u, v, z )dx . 
Jr 

For the stress function a € C 2 (M) we assume the hypotheses 

a'{u)>0, (ai) 
cr(it) = ±|w| p + o(|u| p ) , as u — > ±00, for some p > 1 , (a2) 
and show the following stability and convergence results: 

Theorem 4.2. Zei T > be fixed and let (u, v, z) be any smooth solution 
of (|4.5|) with in particular u bounded in L°° . Let (u, v, z) be a weak, entropy 
solution of (|4.ip such that ^(0) < +00 and 

F e (it, v, z \u, v, z) — > 0, as a; — >• ±00 . 

7/cr verifies ( jail , fla^j ), i/ien 

<C(^(0)+e 4 ), te[0,T), 

for a given positive constant C depending only on T, the properties of a and 
the limit functions u, v and their derivatives. Moreover, if ^(O) — > as 
e -10, then 

sup W(t) -» 0, as e 10. 
te[o,T] 

The proof employs the relative entropy inequality (|4.6p and proceeds fol- 
lowing Theorems 13.31 and 12. 7\ the details are omitted here. 
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